Specific heat of the Gaussian random bond Ising model on a square lattice 
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The free energy and the specific heat of the two-dimensional Gaussian random bond Ising model 
on a square lattice are found with high accuracy using graph expansion and analysis of continuation 
of the high-temperature series. At low temperatures the specific heat follows the power law C oc 
T'^/^, the result confirmed independently by counting many-body states in finite size samples. The 
interpretation of the free energy in terms of independent droplet excitations gives the density of the 
two-level states, that follows a novel power law p(e) oc ^/e at low energies. 
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I. INTRODUCTION 

The disordered two-dimensional systems are ubiqui- 
tous in the modern experiment, reflecting the current 
technological state of growing films. Put simplistically, 
there are two dimensions to disorder: the strength and 
the frustration it induces in the system, both varying 
with the chemical composition. Accordingly, there are 
several types of order at zero temperature: a homoge- 
neous order, when the frustration is small, and an alter- 
native order, the random glass state, frequently observed 
at low temperaturesi*^. In magnetic spin systems with a 
phase transition in a low-temperature ordered phase, like 
the ferromagnet, a weak disorder suppresses the phase 
transition temperature. In strongly disordered systems 
the phase transition in the spin glass state, that normally 
takes place at a finite freezing temperature, Tf, can be 
suppressed by varying the chemical composition, as well 
as by decreasing the film thickness p'"— all the way down 
to zero temperature. At low temperatures the thermody- 
namic properties of the spin glass are determined by its 
low-energy spectrum, described qualitatively by the pic- 
ture of the droplet excitations. In particular, this spec- 
trum determines the lower critical dimension. In this 
paper, for the two-dimensional Gaussian random bond 
Ising model on the square lattice, we find the free energy 
and the specific heat in a wide range of temperatures and 
show, that the density of states vanishes at low energies, 
unlike the density of states in the spin glass. 

The random bond Ising model, describing disordered 
spin systems, is the short-range version of the Edwards- 
Anderson modelii Usually, such models are divided into 
several classes, with continuous and discontinuous, as 
well as symmetric and asymmetric, distributions of the 
bond interactions. For a continuous symmetric disorder, 
of which the Gaussian distribution is a typical example, 
the random bond Ising model develops no order in two 
dimensions'^— with Tf equal to zero. 

In the homogeneous spin systems, where the discrete 
symmetry [Z2 for the Ising system] is breaking down 
spontaneously, the ground state [the pair] is unique, 
and all excitations have a finite energy, the gap. Non- 
frustrating disorder induces tails of the density of states 



inside the gap. In two dimensions these excitations are 
the domain walls, encircling the connected area of the 
flipped spins. The longer the domain wall is, the more ex- 
citation energy it brings in. Similarly, the ground state of 
the two-dimensional Gaussian random bond Ising model 
is probably also uniquei^, whereas the low-energy exci- 
tations are dropletsi^"— , representing connected clusters 
of flipped spins with fractal like boundaries. Since the 
two dimensions lie below the lower critical dimension of 
the spin glass phase^i the energy of the droplets, albeit 
distributed randomly, scales with the droplet size L as 
e (X where the exponent 9 « 0.29 has been found in 

multiple studiesi^ii^iii The boundary of the droplet (the 
domain wall) treads a long fractal path before its total en- 
ergy, locally positive and negative due to the frustration, 
will be fine-tuned into some positive value close to zero. 
From the thermodynamic point of view these droplet ex- 
citations represent two-level, fiipped/non-fiipped, states, 
presumably non-interactingJ^ They determine the free 
energy and the specific heat at low temperatures. The 
extended low-energy excitations can be studied by meth- 
ods of classical statistics, despite that at low tempera- 
tures the quantum fluctuations dominate on small scales. 
In the work [ll| the density of states for a continu- 
ous and symmetric distribution of the bond interactions 
was found to be flnite at zero energy and growing lin- 
early with the excitation energy. Accordingly, the spe- 
cific heat was found to behave linearly, C oc T, at low 
temper ature a ^ ^ 1 ^ ^ . 

In this paper we study the two-dimensional Gaussian 
random bond Ising model on the square lattice at low 
temperatures. The usual method of choice for study 
such models is the Monte-Carlo simulations. However, 
the Monte Carlo method at low temperatures suffers 
markedly from slowing downi^ Owing to a lack of the in- 
tervening phase transition from the paramagnetic phase 
in the spin-glass state, the high-temperature series, de- 
fined and valid in the paramagnetic phase at high tem- 
peratures, can be continued all the way down towards 
the zero temperature, provided a long enough series is 
known. The high-temperature continuation is convenient 
in that it gives a properly disorder-averaged observable 
directly in the thermodynamic limit. We find the free 



energy and the specific heat. The specific heat vanishes 
at low temperature, in accordance with the third law of 
thermodynamics, and is found to follow the power law 
at low temperatures: C(T) (x T", with a sa 3/2. Inde- 
pendently, we check this power law by directly counting 
all the many-body states of the Gaussian random bond 
Ising model on the square lattice, confined onto small 
6x6 samples. The 3/2-law have been also observed in the 
experimentSfii^ though, the origin might be different and 
due to the spin-waves. Interpreting this result in terms 
of independent two-level droplet states, we find, that the 
density of states follows the power law, p(e) oc y/e, at low 
energy e. 

This paper is organized as follows. In section II the 
Gaussian random bond Ising model on the square lat- 
tice is described and the graph expansion method is ex- 
plained. In section III the result for the average free 
energy and the specific heat are found, using the contin- 
uation of the high-temperature series as well as, indepen- 
dently, by studying small samples of the same model. In 
the approximation of independent two-level droplets the 
density of states is found at low energies. 

II. MODEL AND GRAPH EXPANSION 

The two-dimensional Gaussian random bond Ising 
model on the square lattice describes an anisotropic mag- 
netic system, consisting of Ising spins, characterized by 
binary value Sx = ±1, assigned to each site x on the 
square lattice. The two spins across the bond (xy), a 
pair of the nearest-neighbor sites x and y, arc coupled 
by an anisotropic exchange interaction. The total energy 
of this system reads 

H = - JxySxSy (1) 

(xy) 

The nearest-neighbor interaction Jxy is a quenched ran- 
dom variable, which changes from one bond to another, 
and is distributed independently, randomly and symmet- 
rically around zero: 

The average strength of the bond interaction can be nor- 
malized as: J = 1, without any loss of generality. The 
randomness of the bond interactions, Eq.([2]), introduces 
a certain amount of frustration to the spin order, and the 
proper order parameter of this system is expected to be 
of the Edwards- Anderson type: 

QEAiT) = ( (,Sx)|),7, (3) 

where (...)t is the Gibbs canonical statistics thermal 
average with the Hamiltonian, Eq.([T]), at temperature 
T, and (...),/ is the quenched disorder average with the 
distribution Eq.([2]). However, for the two-dimensional 
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FIG. 1. The diagram classes, that contribute to the free en- 
ergy. 



random bond Ising model with a continuous distribu- 
tion of the bond interactions, like the Gaussian dis- 
tribution, there is no Edwards- Anderson order at any 
temperature^'"— qEAiT) = 0. This property is due to 
the presence of the soft modes, representing, for the Ising 
systems, a flipped, with respect to the ground state, clus- 
ters of connected spins. Also, some sites may be con- 
nected to the outside by couplings that all turn out to 
be small, or the "effective" field from the outside spins 
may turn out to be small too. This should be contrasted 
to the homogeneous Ising model on the square lattice, 
where any excitation with respect to the ground state 
carries an energy, greater than the gap A = 8 J, and 
the specific heat follows the Arrhenius activation law: 
C(T) K (/3A)2 exp(-^A), at low temperatures T =1/13. 

We calculate the free energy of the Gaussian random 
bond Ising model on the square lattice, Eq.([T]), using 
the continuation of the high-temperature series. In the 
work [2l] | the high-temperature series for this model was 
calculated, using the linkcd-clustcr expansion method. In 
this method the thermal fluctuations, contributing to the 
free energy in the given order of the inverse temperature, 
/3, are sorted out by their footprint on the lattice, i.e. the 
minimal cluster in which a given fluctuation process may 
occur. It further turns out that many apparently differ- 
ent clusters give identical high-temperature scries con- 
tributions, and therefore can be combined into broader 
classes, represented by a single graph, embeddable in the 
lattice. For that reason we will call this expansion a 
graph expansion. The main novel approach of this paper 
is to re-sum the high-temperature series before summing 
up the graphs. Namely, for each given graph we aim to 
account for all possible fluctuations in all orders of /3, i.e. 
we try to flnd the accurate free energy for each graph in 
a wide interval of temperatures {3. In this way we or- 
der the thermal fluctuations by widening the footprints 
rather than by growing the perturbation order. 

The average free energy per one site of the square lat- 
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tice is given by the graph expansioni22i2^ 



(4) 



fcGa 



where a, 6 numerate the connected simple graphs, that 
are embeddable into the square lattice, ga is the number 
of the embedding of the graph a into the square lattice. 
It comprises the degeneracy under the action of the point 
group D4 of the square lattice on the embedding cluster 
and the degeneracy under the action of the bending parts 
of the embedding cluster. If 6 is a subgraph of graph a, 
then gab is the number of the embedding of the graph b 
into any particular variant of the embedding of the graph 
a into the square lattice. Fa{P) is the average free en- 
ergy of the Gaussian random bond Ising model, limited 
to some embedding of the graph a into the square lat- 
tice. Both gab and Fa{l3) are independent of a particular 
variant of the embedding. The cumulant of the graphs 
a contribution to the free energy is defined recursively 



as 



22.23 



SFaiP) = FaiP) ~Y.9ab5Fb{P). 



(5) 



&Ga 



Note, that Eq.Q follows from the above expression when 
a covers a large part of the square lattice and, therefore, 
asymptotically 5 Fa = 0. 

Let us redefine, shifting by a constant, the bond energy 
of the Ising model as: 



H 



(xy> 



(JxySxSy -Tlog [2cosh(/3Jxy)]) . 



(6) 



Wc observe, that for any graph, that has at least one 
dangling edge, ending in a vertex of incidence valence 
one, the average cumulant of the free energy, evaluated 
using Eq.®, is exactly zero. Eliminating such dangling 
bonds one by one, we end up with a backbone graph, 
that consists of vertices of incidence valence two, three 
and four, all embeddable into the square lattice. Graphs, 
that contribute the cumulants to the average free energy 
of the Gaussian random bond Ising model, Eq.Q, are 
conveniently classified, combined, into diagram classes, 
with the simplest diagrams schematically shown in FiglT] 
These diagrams consist of vertices of incidence valence 
three and four connected by lines that represent vertices 
of incidence valence two and adjacent edges. Obviously, 
the average free energy may depend only on the total 
length of the line, and not on its particular embedding 
route. We also observe, that a cumulant contribution 
from the graphs, that represent one-line reducible dia- 
grams, also vanish. As well as the cumulant contribution 
from the graphs, that represent diagrams, that can be 
disconnected by cutting out a four-fold vertex, the so- 
called articulate graphs. Finally, the sum in the graph 
expansion, Eq.(j4]), span only over graphs from the dia- 
gram classes, shown in FiglT] and the likes. 

Any graph has a distinct property, namely the weight, 
defined as the number of the edges of this graph. Wc 



1.4 



1.38 



1.36 



1.34 



1.32 



1.3 





1 
































beta 







1.5 



2.5 



FIG. 2. The average free energy, —Fw{l3), in units of J, 
vs the inverse temperature, J/3, for four sets of graphs with 
weights up to W=10,ll,12,13 of red, green, blue and pink 
lines correspondingly. Fva^P) is fitted in the interval 0.6 < 
J/3 < 2.2, the black line. (Color online.) 



enumerate all backbone graphs without dangling edges 
up to the weight 14 and calculate their lattice degeneracy 
constants g^, as well as their sub- graph number table gab- 



Ill. SPECIFIC HEAT AND DENSITY OF 
STATES 

We study the low temperature phase of the Gaussian 
random bond Ising model, using a method of continu- 
ation of the high-temperature series, since there is no 
phase transition for the two-dimensional random bond 
Ising model. Summing up in Eq.(j4|) the contributions 
of all the graphs up to weight thirteen, falling into dia- 
gram classes like in Figll] we find the high-temperature 
series for the average free energy per site of the Gaussian 
random bond Ising model on the square lattice: 



F(/3) = - log(2)T - /3 + i/33 - + g/3^ - ^/39 
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(7) 



where the energy unit J = 1 can be reinstated by dimen- 
sion counting. This series is impractical to analyze at 
low temperatures as it is too short and irregular. Instead 
we evaluate the free energies for all the graphs, that have 
weight less than thirteen, more precisely, namely up to 
the 127th order of the inverse temperature f3. This longer 
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high-tcmperature series is re-summed, using the Borel 
method, and then the Pade approximation to the inte- 
grand. Calculating the Borel integral numerically gives 
us a useful approximation of the average free energy de- 
pendence on the inverse temperature. In the interval 
< j3J < 3, thus calculated average free energy for 
any particular set of graphs turn out to be almost ex- 
act. We notice, that rearranging graphs by weight up to 
W results in a partial sum of graphs: Fw{(3)- The high- 
temperature series of Fvi/(/3) in (3 coincides with that of 
Eq.d?]) up to the power 2W — 1, and deviates in higher 
orders. At low temperatures, Fn/(/3) start to increas- 
ingly overshoot in opposite directions, for odd and even 
weights W. And we notice, that by attenuating the con- 
tribution of the subset of graphs of the largest weight W 
by a factor R results in the average free energy: Fwm{13), 
that continues to a somewhat lower temperature. 

We evaluate the average free energy, Fw,r{I3), for four 
sets of graphs of the maximum weight VF = 10,11,12,13, 
all in the interval of the attenuation factor: 1/4 < i? < 
1/2. For each value of W and R an approximation of 
the free energy in the interval 1/2 < /3J < 3 is found, 
using the high-tcmperature series continuation up to the 
127th order. It turns out to be already converged enough 
to be almost independent on the increasing of the max- 
imum 127th order further. Then, the attenuation fac- 
tor is adjusted to give the free energy approximations, 
Pw,B.{P), to the lowest possible temperatures. We find 
R = 0.35, R = 0.40, R = 0.48 and R = 0.36 for the max- 
imum weight, W = 10,11,12,13, correspondingly. The 
resulting approximations, Fw,R{f3), are shown in Figl^l 
For odd weights, W = 11, 13, there is apparently a high 
degree of convergence, whereas curves for even weights, 
W = 10, 12, are less convergent. The anomaly of the 
curve = 12 is coming from one special "window" 
graph. For higher weights the contribution of one graph 
is smoothcned by the others. We assume, that F\y=i3{P) 
is the most accurate approximation of the true average 
free energy per site of the Gaussian random bond Ising 
model on the square lattice in the interval < /3J < 3. 

We fit Fi3{l3) in the interval 0.6 < /3 J < 2.2, using the 
following extrapolation formula: 

F{T)^~{{AT)-' + \E\'')^^\ (8) 

It extrapolates between the ground state energy: F{T) = 
at zero temperature, and the high-temperature be- 
haviour: F = — log(2)T, corresponding to a bunch of 
independent spins. We find, by best fitting 161 points 
of Fi3(/3) with the overall accuracy of 0.02%, that i) 
A ~ 0.73(1), which is indeed close to log(2), ii) 7 ~ 
2.47(6) and iii) the average ground state energy per site 
E = —1.309(3) J, in agreement with what has been found 
before»i^ It is necessarily negative and somewhat higher 
in comparison with E = —2 J of the homogeneous Ising 
model. This effect is due to the frustration and can be 
understood as follows. All bonds on the square lattice can 
be divided into two types, relaxed and frustrated ones, 
according to the value of sign( J^xy)0'xO'y), where Cx is 
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FIG. 3. The density of states in units 1/J vs the droplet 
energy in units J, as found for the Gaussian random bond 
Ising model on the square lattice (red line), vanishing at zero 
energy, and the Ansatz for a spin glass,— which is finite at 
zero energy (blue line). (Color online) 

the ground state configuration of the spins. We presume, 
that the frustrated bonds have relatively smaller interac- 
tions I J(xy) 1 1 whereas the relaxed bonds have relatively 
larger interactions |t/{xy)| in the ground state. Also, 7 is 
the exponent in the power law behavior of the specific 
heat: 

C(r) = (7-1)A (^— j , (9) 

at low temperatures. Approximately, the exponent 7 is 
close to 5/2 and the entropy A is close to log(2). 

In statistical physics the free energy, Eq.®, and the 
specific heat, Eq.®, at low temperatures are both re- 
lated to [i.e. can be derived from] the many-body density 
of energy levels, that follows the law: 

^4„,{„.-.(^)"j, ao) 

for small positive £ — Sq, where £q is the total ground 
state energy, TV is the total number of sites, Z is the 
proper normalization factor, Q is some characteristic en- 
ergy scale of the order of J, and 

7-1 

/3 = 1— . (11) 

In order to verify the validity of this relationship between 
the two exponents, we find all the many-body energy lev- 
els for a small cluster of size 6x6 of the two-dimensional 
Gaussian random bond Ising model on the square lattice 
with periodic boundary conditions and average relative 
energies, £i—£o, for the i = 1..200 lowest many-body ex- 
cited states over 200 or so samples of the Gaussian bond 
disorder. In this way we find the average many-body 
density of states, (logdA/" /d£), and fit it with the Ansatz, 
Ea. ([TU)) . The exponent /3 0.58, given by the best fit, 
turns out to be consistent with the exponent 7 ss 2.47, 
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while representing the independent low-energy approach 
to determine the density of states of the frustrated spin 
system. 

Also, the free energy, Eq.®, and the specific heat, 
Eq.(|9]), can be derived on the microscopic level from 
the model of a non-interacting dilute system of two-level 
states, the droplets. An excited state in a system of Ising 
spins is specified by a set of connected clusters of flipped 
spins with respect to the ground state. Mutually dis- 
connected clusters contribute independently to the to- 
tal excitation energy. The excitation energy of a flipped 
cluster is positive, and is the energy of the domain wall 
boundary. The domain wall energy is the sum of the 
energies of the local bonds, positive and negative alike. 
With a small probability, we can find a specific fractal 
like loop with many negative bonds across it and, if the 
spins inside are overturned, the domain wall energy will 
be atypically small, although positive. Such loops are 
rare in the space of random walks, but in real space are 
dense, and the two-level states are also dense, but pre- 
sumably independent. The free energy of a system of 
independent two-level states reads^ 

F{T) = E-tJ^ p(e)log(l + e-^/^)de, (12) 

where p{e) is the density of states. Working out the 
relationship between the free energy, Eq.®, and the 
two- level states model, Eg. p^ . we find the low-energy 
asymptote for the density of the droplet states: 

^^'^ = 7|i?|7-i(l-2i-7)r(7-l)C(7)' ^^^^ 



Approximately, p(e) oc ^/e at a small droplet energy e. 
The density of states, Eq. (fT3|) . is shown in FiglU as well 
as the Ansatz for the density of states of the spin glass.— 
One can easily see why, provided an insufficient accu- 
racy at low energy, or equivalently low temperature, the 
two relationships can be mistaken. The result, Eq. (|13p . 
shows clearly the difference between the Gaussian ran- 
dom bond Ising model on the square lattice and i) the 
same model defined in one dimension, where the density 
of states p{e) — )■ const at e — >■ and ii) other systems, 
that develop the spin glass phase at low temperatures. 



IV. DISCUSSION 

In conclusion, the average free energy and the specific 
heat has been found with high accuracy in a wide range of 
temperatures for the two-dimensional Gaussian random 
bond model on the square lattice. At low temperatures 
both show a power law behaviour, with C(T) cx T^/^. 
We have found evidence, that the density of states of the 
low-energy droplet excitations vanishes at low energies, 
in remarkable contrast to the finite density of states in 
the spin glass phase and the finite density of states in 
the same model in one dimension. However, our method 
does not reveal the entire spectrum, or other properties 
of the low-energy droplet excitations. 
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